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Abstract 

We report a theoretical model for a thermal transistor in dielectric four-terminal nanostructures 
based on mesoscopic ballistic phonon transport, in which a steady thermal flow condition of system 
is obtained to set up the temperature field effect of gate. In the environment, thermal flow shows 
the transisting behaviors at low temperatures: saturation, asymmetry, and rectification. The 
phenomena can be explained reasonably by the nonlinear variation of the temperature dependence 
of propagating phonon modes in terminals. The results suggest the possibility of the novel nano- 
thermal transistor fabrication. 
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Rapid advances in nanometer-sized techniques have made possible the miniaturization 
and integration of electronic devices. Usually these device properties cannot be described 
by classical transport theory, since their characteristic sizes are small in comparison with 
the elastic mean free path between scattering events, particularly at low temperature. In 
this situation, the wave nature of the electrons needs to be taken explicitly into account, 
and electrons do not propagate diffusively instead ballistically. Based on the ballistic and 
phase coherent electron transport, a variety of interesting devices, e.g., the quantum stub 
transistor, the nanotube transistor, and the multi-terminal junctions, have been proposed^. 

As the counterpart of electron transport, mesoscopic phonon transport has been paid 
much attention recently. The new phenomenon of the universal quantum of thermal con- 
ductance at low temperature was predicted theoretically in mesoscopic dielectric systems 
and verified experimentally [3.]. Some investigations of thermal transport properties have 
been done sequentially at the variety of nanostructuresP, O, 0, 3, 0, [o | - However the ther- 
mal transistor actions in mesoscopic dielectric systems, based on quantum interference effect, 
have to our knowledge not been studied either theoretically or experimentally. It should have 
different phenomena and mechanism from that of the work based on classical theory 10J . 

In this paper, we demonstrate a theoretical model of a thermal transistor in dielectric 
four-terminal nanostructures (DFNSs) based on the ballistic nature of phonon transport at 
low temperatures, in which a steady thermal flow condition of system is obtained to set up 
the temperature field effect of gate. Therefore the DFNSs can work as a ballistic thermal 
transistor, with saturation of thermal flow from source to drain for the gate temperature at 
large temperature difference between the source and drain. For the asymmetric structure, 
thermal asymmetry and rectification are showed under changes in the temperature differ- 
ence. A reasonable explanation of the phenomena is given by the nonlinear variation of the 
temperature dependence of propagating phonon modes in terminals. The possibility of the 
experimental observation of these novel phenomena is also discussed. 

A two-dimensional DFNS is illustrated in FigfjTa). Four parallel terminal wires of width 
cii are directly coupled to a central ballistic window region V with width b and height L 
(scattering region). Other ends of four-terminal wires are connected to thermal reservoirs 
at equilibrium with temperatures T iy i = 1,2,3,4, respectively. Supposed that there is no 
phonon scattering inside the terminal wires and the central window, and a perfect connection 
between terminal wires and reservoirs, phonon scattering is solely decided by the geometrical 
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features and happens only at interfaces between the terminal wires and the central window. 

n 

The net thermal flow Qi in terminal wire i is expressed as[6|] 

Qi = ^2 Qij = ^ ^ f ^?uj[ni(w,Ti) -nj(w,Tj)]Tji inm (w) , (1) 

where rii(uj,Ti) = [exp(hwfkBTi) — is the Bose-Einstein distribution function of the 
phonons with temperature 7$ in the z-th reservoir, is the Boltzmann constant, and h is 
Planck's constant. tOi m = z ^ is the cutoff frequency of mode m in terminal wire i, and v 
(= 5000 m • sec -1 ) is the group velocity. Tj ijnm (uj) is the phonon transmission coefficient 
that an incident phonon with energy hui from terminal i at phonon mode m is scattered to 
terminal j at mode n, with the property: Tji >nm (io) = Tij tTim (u). For the present DFNSs, 
transmission coefficient can be obtained by the mode matching method, the same approach 
as Ref. 0, Q|. Thermal flow is decided by two thermal reservoirs % and j, and has the 
property: Q id = -Q yi . 

In the paper, we assume that the dielectric four-terminal system has reached the steady 
state. The steady thermal flow across central window from terminal 1 to terminal 4 has 
the definition: Q\ = —Q^, where the plus/minus sign represents separately the thermal flow 
into/out of system. Meanwhile temperatures of four thermal reservoirs are assumed to have 
the sequence: T\ ^ (T 2 ,T 3 ) ^ T 4 . In terms of Eq.([T]) and the steady flow definition, we 
readily obtain the steady thermal flow condition of system: 

Q12 + Ql3 + Qi2 + Q 43 = . (2) 

If adding the (Q23 + Q32) into Eq. (J5J), we also get Q2 = —Q3, the steady flow between 
terminals 2 and 3. Obviously, the temperatures of four thermal reservoirs, which are satisfied 
with Eq.(j2j), can make the DFNS at steady state. Note that the summation of four terms in 
the left hand side of Eq.(j2]) is just the thermal flow Q L along the vertical direction of central 
window (region V), which equals zero in terms of Eq.(j2]), (see in FigH^a)). This means that 
the thermal flow across terminals 1 and 4 isolates from that across terminals 2 and 3 at 
steady state of system. With the help of Eq.flT]), equation ((2J) transforms to 

[%C*i) + QlsiTx)} + [QUT±) + QUT A )] = [Q S 21 (T 2 ) + Q* 4 (T 2 )] + [Q S 31 (T 3 ) + Q S U (T 3 )] , 

(3) 

where Q|-(T») is the thermal flow that is produced by single thermal reservoir with 
temperature Tj and that flows directly from terminal i into terminal j; Qf-(Tj) = 
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E m /^ m ^hhJrii(uj,Ti)Ti jtrn (uj), T ijim (uj) = Y^n^i^ ~ u jn )T ij>mn (u). In this way, the steady 
flow condition of system of Eq.Q can be expressed by four single thermal reservoirs. We 
take reservoir 1 for the source with T s = Ti, reservoir 4 for the drain with T d = T4, and 
reservoirs 2 and 3 for the gates with T g = T2 = T3 respectively. The simplified form of 
Eq.(j3]) is written as 

+ Ei{T d ) = E s 2 {T g ) + EI(T g ) , (4) 

where E}(T a ) = Q S 12 (T S ) + Q S 13 (T S ) , E°(T d ) = Ql 2 (T d ) +Q a 43 (T d ), E s 2 {T g ) = Q s 21 {T g ) + Q s 2A {T g ), 
and E^iTg) = Qli(T g ) + Q s M {T g ). Equation (jlj) sets up the source-gate-drain (SGD) relation, 
in which the gates connect the source and drain respectively and work as controlling ends. 
For a fixed gate temperature T g , a pair of corresponding temperatures T s and T d can be 
obtained from Eq.(j4]). Each set of temperatures including the T g and a pair of T s and T d can 
make sure the system at steady state and be used to calculate the steady thermal flow Q sd 
from source to drain in terms of Eq.flTJ); Q sd = Qi in terminal 1 and Q sd = —Q4 in terminal 
4. As stated above, we can get all values of the steady thermal flow Q sd corresponding 
to different temperature pairs of T s and T d that set up the temperature difference between 
source and drain, i.e., T sd = T s — T d , for the fixed T g in the DFNSs. The schematic diagrams 
of working procedure for Q sd is shown in Fig{T](b) — (d). Furthermore, as Q L = 0, the 
influence of the gate temperature T g on the steady thermal flow Q sd should be attributed 
to the temperature field effect (TFE). Note that all temperatures in the paper always mean 
those of thermal reservoirs. 

For the symmetric DFNS, the inset(a) of Fig. [2] shows the calculated thermal flow Q sd 
as a function of the temperature difference T sd = T s — T d for five different values of T g . As 
seen in the figure, Q sd first increases, and then tends to saturation with T sd increasing for 
the five T g s. In the following, we will study the phenomenon based on mesoscopic ballistic 
transport theory. Figure [3] illustrates the cutoff frequencies of the discrete transverse phonon 
modes (threshold energy) respectively in terminal 1, scattering region V, and terminal 4, 
where the superscript p and n represent separately the positive T sd at T s > T d and the 
negative T sd at T s < T d . For the symmetric DFNS with 20nm, 
as shown in Figj3](a), the mode spacing Aj = u^+i) — ^im = ^7 in terminal 1 equals 
that of terminal 4. At initial state with = Tj = T g , the system is at equilibrium 
(Tf d = 0, Q sd = 0). Then the enhance of Tf relative to the fixed T g excites the additional 
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phonon modes AmJ in terminal 1 and makes energy into terminals 2 and 3 from all modes 
{m^ + Amj) in terminal 1, in terms of 7TO ^'" 1 ^ hs |Tf — T g \/h, where the superscript eq 
represents the equilibrium state. Since QfATi) in Eq.flBJ) is monotonously increasing function 
of temperature Tjjf], Q], E{(T S ) in Eq.Q increases correspondingly Meanwhile, to balance 
the both sides in Eq.pj, El(Td) must have equal amount decreased because i£J(T s ) and 
E^(T g ) in Eq.(j4]) keep constant. This means that the modes m e ^ in terminal 4 have AmJ 
decrease and the corresponding temperature T d drops below T g as shown in Figj3]^a). Since 
Amf excited by Tf at Tf > T g are higher energy modes, the decreased energy of AmJ 
which equals Amf is not sufficient to cancel out the additional energy produced by Tf . 
Therefore AmJ in terminal 4 must be more than Arrij in terminal 1, i.e., AmJ > Am{, 
in order to make the increased energy from (m^ 9 + AmJ) in terminal 1 and the decreased 
energy from Am^ in terminal 4 equal. This means that there is an unequal amount change 
of Amj and AmJ for the equal amount of E((T S ) increase and El(T^) decrease in Eq.([4"l). 
With Tf d increasing further, as Amf ^ \Tf — T g \, the ratio of AmJ to AmJ is the nonlinear 
monotonously decreasing; therefore Tf shows the nonlinear increasing behavior, meanwhile 
Tj has the nonlinear decreasing behavior as shown in FigH^a), in order to keep the system 
at steady state. Obviously the modes (AmJ + AmJ) corresponding to (Tf — T d ) become 
more with Tf d increasing, but the ratio of Am^ to (AmJ + Am^) goes down. This means 
that Tf has less contribution than T d for Tf d increasing in terms of Amf ^ \Tf — T g \. On 
the other hand, as the propagating channel between the source and drain consists of the 
propagating phonon modes (m^ 9 + Awij) in terminal 1 and (m^ 9 — AmJ) in terminal 4, Q s d 
increases with Tf d ; but the increment rate of thermal conductance Q s d/T^ d reduces with 
Tf d increasing, i.e., Q s d/T^ d > and -^(Qsd/T^) < 0, (see in Figfjjb)). Consequently, 

sd 

under the influence of T g , the nonlinear variation of the modes AmJ and AmJ in the modes 
(Amf + AmJ) that corresponds to Tf d leads to the saturation behavior of thermal flow 
Q s d in the DFNS. Furthermore, as the mode spacings in terminals 1 and 4 are identical 
in this symmetric structure, Q s d presents symmetry in the reversal of Tf d , i.e., T d > Tf . 
Additionally, Q2{= —Q3) is two orders of magnitude lower than Q s d- 

In order to investigate the asymmetric effect, we change the DFNS from the symmetric 
structure to the asymmetric with = 20nm, 04 = 40nm. As shown in FigfSJ there 

is asymmetry in the Q s d— T s d behavior with thermal flow saturation for the positive and 
negative T s d for four values of T g . (i) Thermal saturation discussion: for the asymmetric 
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DFNS, the mode spacing Ai = ^ in terminal 1 is twice A 4 in terminal 4; there is a 
higher density of transverse phonon modes in terminal 4 as shown in Figj3](b). When Tf 
increases over T g , thermal energy input from the modes (mJ'+AmJ) in terminal 1 is canceled 
out by the decreased modes AmJ along with T d lowing, to satisfy the SGD relation, (see 
in Figj3]^b)and FigJUc)). However, because there is higher density of phonon modes in 
terminal 4 than that in the symmetric situation, at initial stage, the small change of T d 
could make sufficient modes AroJ decreased to cancel out the energy input from Tf enhance; 
the ratio of AmJ to (AmJ + AmJ) drops less with T v ad increasing, in comparison with the 
symmetric situation. Therefore Q s djT v sd increases with T v sd increasing, i.e., Q s d/T^ d > and 
dT p ~(Qsd,/Ts d ) > 0, (see in Figfjjcl)). With Tf increasing further, more AmJ with higher 

sd 

energy in terminal 1 are excited and more energy is input into system; therefore AmJ in 
terminal 4 reduce more and T d correspondingly goes much down to balance the both sides 
in Eq.(j4l). The ratio of AroJ to (AmJ + AmJ) decreases remarkably with Tf d , similar to the 
symmetric situation. Thus Q s d/T^ d changes from increase to decrease with Tf d , as shown 
in Figged); -^(Qsd/T^) > at beginning, and then -Ar{Q s d/T^ d ) < after T sd passing 

sd sd 

a critical value. This means that Q s d — Tf d curve increases monotonously with the concave 
shape at beginning, and then keeps going with the convex shape after passing a turning point, 
and finally tends to saturation, (see in Figj2]). (ii) Thermal asymmetry and rectification 
discussion: in the reverse case with T d > Tf , due to the higher density of phonon modes in 
terminal 4, T d enhance over T g could easily excite more Am™ in terminal 4 and makes more 
energy into system. This leads that Am™ in terminal 1 must decrease much, due to the 
larger mode spacing A = ^ in terminal 1. Thus Tf corresponding to Am™ decreases more, 
in terms of 7TV ^ 1 ^ &b \T™ — T g \/H, (see in Figj3](b)). The ratio of AmJ to (Am" + Am") 
decreases rapidly with T™ d , so that Q sd tends more quickly to saturation in the reverse case 
as shown in FigfJl Furthermore for the fixed T sd : \T S — T d \ = (Tf — T d ) = (T d — T") = 
constant, the AmJ with higher energy for the positive case (Tf > T d ) is more than the Am" 
for the reverse case (Tj 1 < T d ); therefore we have Tf > T d and T d > Tf , (see in Figj3](b)). 
Thus, in the positive case, Tf and T d locate at higher temperature region to make more 
modes with higher energy participated in thermal transport; for the reverse case, Tf and T d 
is situated at lower temperature region to excite less modes with higher energy for thermal 
transport, (see in FigJD^c) and Figj3]^b)). Consequently the Q s d — T sd characteristics display 
a distinct asymmetry and rectifying behavior in the asymmetric structure as shown in Figj2l 
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In the above discussion, the main effect of T g is to set up a reference temperature which 
makes the nonlinear variation of Ami arid A7714 possible. 

In view of some successful thermal measurements in nanoscale , it is possible 

to realize the relevant experiments for the novel phenomena in the DFNSs. Authors think 
that the steady thermal flow in the ballistic system is a necessary factor to observe these 
thermal transisting behaviors as discussed in the paper. 

In summary, thermal transistor has been designed theoretically in the DFNSs. When sat- 
isfied with the steady thermal flow condition of system, saturation and asymmetry and rec- 
tification of thermal flow are exhibited separately in the symmetric and asymmetric DFNSs. 
The relevant physical mechanism is elucidated by the nonlinear variation of temperature 
dependence of phonon modes in terminals. The mechanism is also helpful to investigation of 
the negative differential thermal resistance and thermal diode design in nanoscale systems. 

This work was supported through grants by the Hong Kong Research Grants Council 
(RGC) and the Hong Kong Baptist University Faculty Research Grant (FRG). 
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FIG. 1: Schematic diagrams for the research system and the working illustration of SGD. Arrow 
lines represent the thermal flows, (a) A four-terminal structure, (b) SGD at equilibrium: T s = 
T g = Tj. (c) SGD at steady state: T s > T g > Tj; a propagating channel (dash lines) is set up 
between the source and drain, (d) Thermal flow profile at steady state. Q\ is the net thermal flow 
in terminal 1, and Q L is the net thermal flow across the level line (dot line) or (arrow line) in the 
vertical direction of region V and always equals zero at steady state of system. 
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FIG. 2: Thermal flow Q sc i vs temperature difference T S( i for the gate temperature T g in the DFNS, 
the same structure as shown in FigQJa). [Main figure] Asymmetric structure with = 40nm. 
[Inset (a)] Symmetric structure with 04 = 20nm. Other parameters: a% = a>2 = 03 = 20nm, 
L = lOOnm, b = 20nm. 
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FIG. 3: Schematic diagrams for the cutoff frequencies of discrete transverse phonon modes or 
threshold energy (straight lines) for terminal /, region V, and terminal IV. Dash dot dot line: 
the gate temperature T g . Dash dot lines: temperatures T s and T^. Arrow lines: the change of 
temperatures relative to T g . The superscript p represents the positive T S( i (T s > T<j); n for the 
negative T s d (T s < Td). AmJ is modes in terminal 1 excited by Tf. (a) Symmetric structure, (b) 
Asymmetric structure for Xf , = T™ d . 
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FIG. 4: At steady state of system, for the fixed gate temperatures T g s: [(a) and (c)] Temperatures T s 
in source and Td in drain vs temperature difference T sc i between the source and drain, respectively 
for the symmetric and asymmetric cases, [(b) and (d)] Thermal conductance Q s d/T s d vs T s d 
respectively for the symmetric and asymmetric cases. 
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